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Abstract 



In this paper we explore the general conditions in order that a 2-dimensional natural Hamiltonian system 
possess a second invariant which is a polynomial in the momenta and is therefore Liouville integrable. 
We examine the possibility that the invariant is preserved by the Hamiltonian flow on a given energy 
hypersurface only (weak integrability) and derive the additional requirement necessary to have conserva- 
tion at arbitrary energy (strong integrability). Using null complex coordinates, we show that the leading 
order coefficient of the polynomial is an arbitrary holomorphic function in the case of weak integrability 
^ ' and a polynomial in the coordinates in the strongly integrable one. We review the results obtained so far 

, with strong invariants up to degree four and provide some new examples of weakly integrable systems 

with linear and quadratic invariants. 
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H ■ 1 Introduction 

In 1983 Hall [1 published a remarkable paper devoted to a theory of configurational invariants of classical 
Hamiltonian systems. The main virtue of the work consisted in an elegant and powerful technique to solve 
the equations for the existence of a second polynomial invariant of arbitrary degree in the momenta for 
2-dimensional Hamiltonian systems. As a result of this approach. Hall was able to get some new examples 
of integrable systems admitting a second invariant of degree four in the momenta. 

Unfortunately, the paper was flawed by a definitely wrong statement and also by many inaccurate ar- 
guments and deductions. In fact. Hall purported to remedy supposed oversights in previous works on the 
search for the second invariant, in particular criticizing the classical account by Whittaker In Hall's 
view, Whittaker's (and all others' since then) treatment provides only sufficient conditions for the existence 
of a second invariant (linear and quadratic in the specific instance), whereas, due to overlooking the link 
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established on phase-space variables by energy conservation, it was not able to find all possible solutions. 
Actually, as it is, this statement is wrong: in looking for a strong second invariant, namely a phase-space 
function which commutes with the Hamiltonian function, Whittaker's approach is indeed correct and leads 
to necessary and sufficient conditions for its existence. This point was already stressed by Sarlet et al. [3J in 
their criticism to Hall's paper. 

Moreover, Hall's discussion contained a confusion between the concept of configurational (or weak) in- 
variant, as a function which exactly commutes with the Hamiltonian only on a subset (possibly one) of the 
energy hypersurfaces (see again [3] and, e.g. [1]), and the notion of what we may call formal integral as it 
emerges in the analysis of regular portions of the phase space of generic non-integrable systems (see e.g. [5]). 
In particular, the approximate invariants obtained by Hall as a result of his perturbative approach have little 
to do with those that can be obtained by truncating a normal form expansion. 

Nonetheless, in spite of all the above shortcomings, the form in which the problem has been set by Hall and 
the approach followed for its partial solution deserve attention, since they can still be very useful. Already 
Hietarinta f6] , in his account of the direct methods for the search of the second invariant, provides a review of 
all the known systems admitting one or more configurational invariants and works out again the integrability 
conditions found by Hall for the existence of weak invariants up to degree four. However, in Hietarinta's 
review, the two settings of weak and strong integrability are still kept well separated. More recently, in a 
series of works about a unified approach to treat both kinds of invariants [3 [H [9l [14] , the same integrability 
conditions have been obtained, supplemented by the additional constraint imposed by strong integrability. 
This last step is essential for a neat distinction between the two notions of weak and strong integrability. 
In these works, this step is a straightforward consequence of the geometric approach in which the existence 
of the second invariant is addressed by studying the corresponding Killing tensor equations for the Jacobi 
metric. This metric depends on the mechanical energy as a parameter. Therefore, any tensor object on the 
corresponding manifold in general depends on the energy. With due care about formal relations between the 
two approaches, the geometric approach and Hall's approach are equivalent. The essential remark is that, 
to identify the cases of strong integrability, it is sufficient that in the final results concerning the existence of 
the invariant, a subset can be isolated which is independent of the energy parameter. In this framework, in 
[7], quadratic invariants at arbitrary and fixed energy for 2-dimensional Hamiltonian systems were treated 
in a unified way, whereas in [8l [9] the existence of respectively cubic and quartic invariants was discussed 
accompanied with the discovery of some new examples not given in earlier works (see e.g. [TUl [TTl \T2[ IT^ l. 
In [14] the case of Hamiltonian systems with vector potentials is treated, extending previous investigations 

[miiiiniiisiiiniiio]. 

The aim of the present paper is to discuss the techniques for solving the equations for a second invariant 
having a momentum dependence of arbitrary polynomial degree. The analysis is based on a combination of 
both the above-mentioned approaches, with particular attention dedicated to the conformal transformations 
used to simplify the equations. The treatment is in general effective for both classes of invariants, config- 
urational (or weak) and strong. After that, we impose the additional condition needed to isolate the class 
of strongly integrable systems, obtaining the general form of the simplifying family of transformations for 
each degree of the strong invariant looked for. This in turn implies determining the leading order terms in 
the invariant itself. An alternative route to this result in the geometric approach is based on the invariance, 
under conformal transformations, of the conformal part of the Killing tensor. Moreover, a general review of 
the results obtained so far is provided. 

The plan of the paper is as follows: in Section 2 we compare the notion of weak and strong invariants 
and, working out in detail the quadratic case, we correct Hall's misunderstanding, providing the constraint 
to be satisfied in order to get strong integrability as a restriction of weak integrability; in Section 3 we recall 
time reparametrization of the null Hamiltonian and complex ( "null" ) coordinates used to simplify the direct 
approach; in Section 4 these tools are exploited to set the general approach to find polynomial invariants of 
arbitrary degree; in Section 5 we recall the main results concerning invariants of degree up to four; in section 
6 we present our conclusions and the prospects for future works. 
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2 Configurational invariants versus strong invariants 



As it is well known (see, e.g. [21 ), to grant the complete integrability of an A^-dimensional Hamiltonian 
system, it is necessary and sufficient to find N independent integrals of motion or invariants, for short. 
In the following, we will limit ourselves to the simplest case of a conservative Hamiltonian system in two 
dimensions {N = 2). Since the Hamiltonian itself 

H^i{pl+p'y) + V{x,y), (1) 

is a conserved function, it is enough to find just a second independent invariant. 

In his search for polynomial invariants. Hall [T] exploits what is usually called the "direct approach" 
of Darboux [22] and Whittaker [2] . As an "educated guess" we may start with the working hypothesis of 
an invariant with a structure analogous to that of the Hamiltonian, namely a second degree polynomial in 
the momentalll Clearly, one could start with a polynomial of arbitrary degree or even with a more general 
expression. In the core part of the paper we will examine the general case. Now, as an introduction aimed 
at giving the proper settings, we work out in detail the quadratic case. We may then look for a phase-space 
function of the form 

l2{Px,Py, X, y) = Apl + Bp^py + Cpl + K, (2) 
which is preserved along the flow given by eq.lJT]), namely 

{h,H} = 0. (3) 

In eq.(l2]). A, B, C and K are each functions of x and y. Actually Hall, probably motivated by his interest in 
accelerator physics, manages to work out the integrability of Hamiltonians including also a vector potential 
and henceforth terms which are linear in the momenta. This, in turn, suggests the inclusion of analogous 
terms in the invariant. However, in the standard case of a Hamiltonian invariant under momentum inversion, 
p —p, the equations for the coefficients of even and odd degree terms decouple, greatly reducing the 
complexity of the system to solve. Therefore, since the main purpose of this paper is to fully clarify the 
issues of weak and strong integrability, which are unaffected by such generalizations, we prefer to limit 
ourselves to the standard case of eq.((T]). The treatment of the vector potential is presented elsewhere [l4] 
where new strongly integrable systems with quadratic invariants are presented. 

In the direct approach, one inserts the functions (ll|2|) into the Poisson bracket ([3]). The resulting poly- 
nomial in the momenta, of third degree in the present instance, must be identically vanishing. In view of 
the independence and arbitrariness of the momentum coordinates, each coefficient of the polynomial must 
vanish, determining in turn the following system of PDEs in the coordinates: 

A, = 0, (4) 
Ay + B, = 0, (5) 
By + = 0, (6) 

Cy - 0, (7) 

K, = 2AV^ + BVy, (8) 
Ky = BV, + 2CVy, (9) 

where, in order to compactify expressions, suffixes denote partial differentiation. As mentioned above, 
Darboux [22] addressed the problem and the equation ensuing from the integrability condition for K, that 
is 

B{Vyy - V^,) + 2(A - C)V,y + ^iByVy - B .,V^) = 0, (lO) 

is known as Darboux's equation. Whittaker [5] also analyzed the problem and, for a complete account of the 
solution of this system, the standard reference is again the review by Hietarinta [B] . 



Since all the work is made in a Hamiltonian context, contrary to the original treatment by Hall, we use canonical phase-space 
coordinates. However, we have tried to stay as close as possible to his treatment as refers to the procedure and notations. 
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Equations involving only the leading order terms {A, B and C) in the invariant are readily solved: 

A^ay'^ + by + c, (11) 
B = —2axy — bx — dy — e, (12) 
C = ax^ +dx + f. (13) 

Exploiting linear transformations of the coordinates, it can be shown [51 [7] that these functions can be 
reduced to four canonical forms that, inserted into Darboux equation pO)) . lead to the complete solution 
of the problem: there exist four fundamental separating coordinate systems (elliptical, polar, parabolic and 
Cartesian) in which the potential and the remaining unknown of the invariant, can be expressed in 
terms of combinations of arbitrary functions of each of the separating coordinates. Actually, in his reference 
to Whittaker's result. Hall mentions only the elliptical solution, but, as shown in [7] (see also [13]), using 
conformal transformations, all the separable cases possess the same structure. Moreover, we remark on 
the additional possibility of complex potentials, which provides four more cases with complex separating 
coordinates, which are listed in [3] and can also be obtained with the techniques in [7]. 

At this point Hall states that ([T], p. 93) "Whittaker's work /.../ was flawed by his failure to recognize 
that terms in and are not independent for the purpose of setting the coefficients of their various 
powers to zero /.../ They are related by conservation of energy. Therefore, Whittaker's constraints were 
over-restrictive; the solutions he found are valid, but others may also exist." On these bases. Hall elaborates 
a generalized approach that, in his view, provides necessary and sufficient conditions for the existence of 
the given invariant. Actually, as we will shortly see. Hall fails in recognizing that the possible extension of 
the family of solutions implies a different status for the additional invariants: they are "fixed energy" or 
configurational invariants. 

Configurational (or "conditional" or even "weak" ) integrals hold for a specified particular value of the 
energy constant and were firstly investigated by Birkhoff 24J and considered by Fomenko ^25], Kozlov [26j 
and others. Weak invariants enjoy weaker properties then their "nobler" cousins, the strong invariants, 
which are constant on every energy hypersurface admitted by the dynamics of the system. To clarify this 
essential point, let us again work with the quadratic case, concretely introducing Hall's argument. This goes 
as follows: for every value E of the energy, the function ([T]) defines the hypersurface 

l{pl+pl) + V{x,y)^E. (14) 

Let us construct the following linear combinations of the squares of the momenta: 

(15) 

A=pI~pI (16) 
and observe that, in view of (|14|) . we can impose the constraint 

S = 2G(x,2/), (17) 
where we have introduced the "Jacobi" potential 

G = E-V. (18) 
The quadratic terms in the standard form of the invariant of eq.Q can now be written as 

1 r 



Apl + Bp.,py + CpI^- \ {A + C)S + (A - C)A + B^Jy? - A2j . (19) 
Exploiting the energy constraint (|17p and redefining coefficients, the quadratic invariant can then be written 



as 



l2{Px,Py, X, y) = \D{pI - pI) + Bp.py + K, (20) 



^The origin of this denomination comes from the close relationship between the picture of a system constrained on the fixed 
energy surface and the geometric picture of a geodesic flow over a Riemannian manifold endowed with a "Jacobi" metric (see, 
e.g. [27]). 
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where 



D = A~C, (21) 

K^K + {A + C)G. (22) 

Now, we may proceed along the same hnes foUowed above. The commutation relations ([3]) of the weak 
invariant (|20p with the Hamiltonian reduce to the system 

Dy + S, = 0, (23) 

D,. - By = 0, (24) 

= BVy + DV., - GBy, (25) 

Ky = BV^~DVy-GB,. (26) 

The integrability condition for K is now the generalized Darboux equation 

B{Vyy - 14,) + 2DV,y + 3{ByVy " B ^V^) - 2{E - V)D.,y = 0, (27) 

where we have explicitly pointed out the presence of the energy parameter. In practice, since the energy of 
the system is in general an arbitrary real number, we can write eq. ()27p in the form 



fiE + /o = 0, (28) 

where /i, fo are functions of the coordinates, both explicitly and through V and its derivatives. 

If we want that the looked for invariant has an identically vanishing Poisson bracket with the Hamilton 
regardless of the value of the energy (in other words that I2 should be a strong invariant), eg. ([28]) must be 
satisfied for every value of the parameter E so that the two equations 

fi - 0, (29) 
fo = 0, (30) 

must separately be satisfied. The first of these equations is simply 

D^y = 0. (31) 

As a consequence of this, the second, by a direct comparison, turns out to coincide with the standard 
Darboux equation (fTU|) . This time, the solution of equation ([3T|) . together with (^3)1 and ([M)) for the leading 
order terms (now D and B), gives 

D = a{y'^ ~ x'^) +by - dx + g, (32) 
B — —2axy — bx — dy — e. (33) 

Recalling ((2T|) and comparing with (fTT|[T3|) . we see that we have arrived at a result completely equivalent to 
that of the standard Darboux- Whittaker approach. In fact, cq. ((30|) now coincides with cq. (fTO|) and, therefore, 
the same possible separable potentials in 2 dimensions can be found also following this alternative route. 

In our general presentation in section 4, we will see how this strategy reveals to be useful also in the 
case of higher degree invariants. While used mostly for its pedagogical role in the present instance, what is 
important to stress is the key role of the integrability condition in the form ([28]): in the general case of a 
polynomial invariant of degree M, it turns out that eq. ([28|) is a polynomial in E of degree M/2 for M even 
and (M + l)/2 for M odd. 

In the case of invariance at fixed energy, the problem admits additional solutions. Now we have that 
ea. ([3T|) is no more necessary and therefore ([^3]) and ([M]) are unconstrained Cauchy-Riemann equations. 
Instead of (|32p and l|33[) . the solution is now given by an arbitrary analytic function: its real and imaginary 
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parts respectively provide the functions D and B. Because of this feature, we will see in the following how 
it is more convenient to work with complex variables. 

Summarizing, the above procedure shows that the strategy for finding weak invariants is well defined and 
points out where it must be constrained to get strong invariants. This objective being the most important 
for applications, it may appear that, going this way, nothing is gained with respect to the usual direct 
approach. However, the procedure instead proves to be very effective in simplifying the system of equations 
resulting from implementing the direct approach. This simplification is actually one of the reasons why 
Hall's work is still useful. At the same time, since the approach based on the constrained invariant ansatz 
(PU|. complemented by a correct use of the general Darboux equation (|28p . provides the same results as the 
standard approach, we get a simple proof of the invalidity of Hall's criticism towards Whittaker. However, 
Hall is right in envisaging additional solutions to the problem: in the example we have just seen, it is 
conceivable to obtain, for particular values of solutions of eg. ([28]) not included in those of[29land[29l In 
[7], we have produced several classes of solutions corresponding to i? = 0. 

The investigation of such weakly integrable systems (WIS) offers several additional issues to study. We 
mention some of them: 

- possible solutions defined for a continuous but finite range of energy values Ei < E < E2; 

- WIS which are actually SIS (strongly integrable systems) with a more general form for the second invariant 
(Examples: the Kepler problem, the Sarlet et al. case [5); 

- the main property of the integrable dynamics on the E surface may help in understanding some of the 
main features in the non- integrable regime (see e.g. [55]). 



3 The null Hamiltonian and time reparametrization 

As already put forward by Hietarinta (see [5], sect. 7. 2), canonical point transformations generated by ana- 
lytical functions preserve the form of a null Hamiltonian. This invariance allows a straightforward way to 
reduce the set of equations to be solved in the search for a polynomial invariant. In the present section we 
show how, in the case of 2-dimensional systems, the above transformations are actually conformal trans- 
formations and are related with the time reparametrization of the dynamics. These tools, exploited in the 
remaining part of the paper, were implicit in Hall's work. 
In general the Hamiltonian itself has the form 

H = T + V{q)=E, (34) 

where T is a quadratic form in the momenta. The independent variable, let us say t, is often but not always 
the time. For any given energy E of the system, to represent the dynamics, we can use the null Hamiltonian 

Ho = H -E, (35) 

provided that we impose the constraint 

Ho = 0. (36) 

For any such zero energy Hamiltonian we can reparametrize the system by introducing a new time variable 
t defined by the relation 

dt = N{p, q)di , (37) 

together with a redefined Hamiltonian 

Ho ^ N{p, q)Ho ^ NT + N{V ~ E) ^ . (38) 

The new Hamiltonian will then give the same equations of motion on the constraint surface Hq — 0. We 
shall use the term lapse function for N{p, q) which defines the independent variable gauge. This usage is 
borrowed from cosmological applications where the lapse gives the rate of physical time change relative to 
coordinate time. The lapse function can be taken as any non-zero function on the phase space. 
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It is simpler to work with complex, or "null" , coordinates, 

z^x + iy, Pz = P ^ iiPx - iPy), (39) 
z^x~iy, Pz^P = h{Px+iPy)- (40) 

The null Hamiltonian can then be written in the form 

i/o = 2pp - G(z, z) = 0, (41) 

where G is the function introduced in (fT8|) . The equations of motion given by ([4T|) are 

§=2p, ,42) 

and corresponding complex conjugates. G and V are always assumed to be real functions. In the following 
we will spare to mention explicitly to the complex conjugates. 

We use a conformal transformation to standardize the frame and coordinate representation of the invari- 
ant. To that end we introduce a new complex coordinate w by means of the transformation 

z = F{w). (44) 

The conformal transformation 

z~^w = X + iY, (45) 
given by the holomorphic function (j44p determines the canonical point transformation 

w^F-^{z), P^F'p, (46) 

so that ([1]) transforms into the new null Hamiltonian 



~ _ 2PP^Giw,w) 
^0 0. (47) 



Let us introduce the "standard" null Hamiltonian 

Hs ^2PP-G{w,w) = 0, (48) 
that implies the use of the new "time" s, such that 

as at 



2P, (50) 
G^- (51) 



Equations of motion ([42H43)l becomes 

dw 

ds 
dP_ 

ds 

At the same time, a conserved quantity stays conserved if transformed between the two gauges (|47|) and 
In terms of real variables, (|48|) is given by 

Hs^\{Pl+Ph-G{X,Y)^Q. (52) 
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4 Invariants of arbitrary degree 



With the choice of polynomial invariants, the time reflection symmetry of Hamiltonian ((1} allows a further 
simplification in the procedure. In fact, it can been proved (see Hietarinta [6], sect. 2. 3 and [22]) that 
the algebra of commuting functions with the Hamiltonian, which is an even function with respect to time 
reflection, has "good" time reflection parity. Therefore, a polynomial second invariant can be either even or 
odd polynomial in the momenta. We can therefore assume a phase-space function of the following form: 

[A// 2] M-2k 

i; = i=0 

where the functions ^(j,M-2fe) (2;, y), not necessarily polynomials, have to be determined. In the first sum- 
mation, with [Af/2] we denote the greatest integer less than M/2, so that if, e.g., M — 1, k takes only the 
value zero. 



4.1 The direct approach in the Cartesian frame 

The simplest procedure (also called the direct approach) is now to compute the Poisson brackets of Im with 
iJ, collect terms with various power of momenta end let vanish their respective coefficients. Using at first 
the form (I53|) in the usual Cartesian frame, we then get a system of partial differential equations of the form: 

(j + l)A(^j+i^k+i)dxV + (fc + 1 - j)A(^j^k+i)dyV = a^r^o-i.fc-i) + dyA(^j^k-i), (54) 

with j = 0, fc and fc = M + 1, M — 1, Af — 3, or 1 and it is implicit that ^(^.f) = if s < or s > i and 
t < OT t > M. In order to simplify formulas, in the present section we reintroduce the standard notation d 
to denote partial differentiation with respect to the variable in the subscript. The set is a complicated 
set of PDEs, in general overdetermined: if M is even we have (M -I- 2)(A/ -I- 4)/4 equations for (A/ + 2)^/4 
unknowns; if M is odd we have {M + 3)^/4 equations for {M + 1){M + 3)/4 unknowns. Actually, in the 
general case, to these figures the potential V enters as an additional unknown. 

4.2 The direct approach in the null frame 

Introducing the complex null frame (I39I40|) . we have the following generic expression of the invariant: 

( M/2 ^ 

hi = 2Rc j C2fep2^ \ , even (55) 

C (M-l)/2 ^ 

hi = 2Re j Yl C2fe+ip''+' \ , M odd, (56) 

where the complex functions C2fe,fc = 0,...Af/2 or C2k+i,k ~ 0,...(Af — l)/2 depend on z,z and, with the 
exclusion of the leading order coefficients Cm, implicitly on E. Their explicit expression in terms of the 
coefficients of the invariant in Cartesian form are 



— -fc 

' 2{a+k) 



X 



■i+k)-i 

i=0 



r=0 



e \/ 2(a + fc)-^ 
e~a + r) \2{a + k)~e~r 



(57) 



A(e,2(a+k)) / , 
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when M is even and 



MM) 




^2(a+fe) + l-C 



(58) 

2{a + k) + l-£ 
e-a + rj \ 2{a + k) + l~£-rj 



^(£,2(a+fe) + l) 



when M is odd. In these expressions a superscript (Af) has been introduced in order to let them apply 
in general: namely, ci^^^ denotes the a — th coefficient of the invariant of degree M. In the following this 



superscript will not be used unless it is strictly needed to prevent confusion. As usual, in ((57)) and (|58 



kj {n-k)\k\ 



(59) 



denotes the binomial coefficient. 

The forms ()55p and (I56p of the invariant do not contain cross terms, namely terms with powers of pp, 
since they have been eliminated exploiting the constraint 

2pp = G{z, z) (60) 

dictated by ([1T|) . A simple check of this statement can be performed with the quadratic invariant with which 
we are familiar from Section 2. Eg. ([57)1 with M = 2 give 

=^22-Ao2+iAi2, (61) 
C^'^ =^00 + (Ao2+A22)G'. (62) 

Posing A22 = A, A12 = B, A02 ^ C, Aqq ^ K and ^ ^ D + iB, ^ K we again find relations ([^11^ 
and see that the expression (j20[) of the invariant is equivalent to 

h {P, P,z,z) = Cf^p" + + K_ (63) 

The use of the energy constraint leads to the maximal reduction in the number of equations ensuing from 
the direct approach. However, in the case of integrability at arbitrary energy, to this reduced number one 
must add the integrability conditions which results in a set of equations generalizing ([29]) . The constraint ([60)1 
can also be used to "homogeneize" the polynomial invariant. In this form, used in the geometric framework 
in coimection with the geodesic flow over a Riemannian manifold, the coefficients of the invariant give rise 
to a symmetric M-rank tensor, known as a Killing tensor. 

The commutation relation of the functions or with the null Hamiltonian (|¥T|) gives the system 
of equations 

d,Ck-i + \{d,Ck+i)G+\{k + \)Ck+id,G^Q, fc = 0,l,...,M, (64) 

where it is implicitly assumed that Cj = for < and for j > M . The set of equation ([64| must be 
supplemented by the closure equations 

d-,CM = (65) 

and 

^{d,{CiG)} = 0. (66) 

Clearly, this last condition must be stated only in the case of odd degree, since for M even, a closure condition 
of the form 

9{Co} = (67) 

is implicit since Co is real by definition. 
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In both cases, using the null complex coordinates, we have M + 2 real equations for M + 1 real unknowns: 
the substantial reduction of the number of independent equations with respect to that in the Cartesian frame 
was already remarked by Hietarinta (see [B], sect. 7. 5). However, he thought that reintroducing the explicit 
dependence on energy would reestablish the original number, whereas, as we will see, the gain in saving 
equation still remains when going to the Cartesian frame. 

4.3 Solving the equations for the Mth-degree invariant 

At present a general solution of system ([M)) is still lacking. In [3 [HI [S] we found the general solutions for 
1 < M < 4. In the present subsection we illustrate the aspect of the procedure which are common to all 
Mth-degree invariant. The first steps are essentially the same as in the original paper by Hall. 
Looking at the system above, we see that eq. (|55|) is readily solved: 



that is Cm is an arbitrary holomorphic function. The first important result we get is therefore that the 
leading order coefficient of a polynomial invariant is given by an arbitrary holomorphic function. This fact 
is already known in the case of homogeneous polynomial invariants (Kolokoltsov [30], Kozlov [26]) and was 
obtained by Birkhoff ([21], chap. 2) for M < 2. We remark that, in agreement with notations in ([57|l 
and (|58p . here we are referring to that is the leading order coefhcient of the polynomial invariant of 

degree M with various values of A/. Therefore, in order to avoid confusion with lower order coefficients in 
a polynomial with a given M and also to conform with the notation in previous works, in the following we 
will denote the leading order coefficient with Sm- 

In order to attack the remaining equations in (I64|) . we may further simplify them by performing a 
coordinate transformation. This is a conformal transformation of the form (|44p . where the generating 
function F{w) is chosen such that 



Cm — Cm{z), 



(68) 




(69) 




dn,CM-2 + ^-^n,G ^ 0, 



(70) 



where, in agreement with (|47p . the new "conformal" potential 



G{w,w) = \F'{w)\'G = [SmSmY'^^'G, 



(71) 



has been introduced. Defining the function 



Z — iVl 

Cm-2{'w,w) ^ Sj^l" Cm -2, 



(72) 



ea.([7n|) becomes: 



dtuCM^2 + MrdwG = 



(73) 



and, defining in an analogous way. 



C2k{w, w) = S^/" C2k, fc = 0, 1, . . . , [A^/2], 



(74) 



the rest of the system becomes 



5ffiC,-2 + -J—d^CjC^) = 0, 



(75) 



with the ranges = 4, 5, . . . , iV — 2 if iV is even and j = 3, 4, . . . , TV — 2 if is odd. 
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4.4 The Kahler potential 



Eg. ([70)1 ■ or better its transformed form ([75]) . plays a special role. In view of the fact that the new potential 
G is still a real function, (|73p can be solved introducing a real function IC{w,w) such that 

G = (76) 

CM-2^-¥dl^lC. (77) 

The first equation says that the Jacobi potential is the Laplacian of the function /C(w,u)). Since, in the 
geometric picture, G is the conformal factor of the metric element, using JC it gives rise to a Hermitian form 
and therefore it is referred to as the Kahler potential. Without loss of generality, in ([2|, sect. IV) it has 
been shown that the Kahler potential can be assumed to be of the form 

K. = E[F{w)F{w) + 25R{A(z(w))}] - (78) 

where A is an arbitrary holomorphic function independent of E and the real prepotential ^I' is such that 

dl^^^F'F'd^,,^ = \F'\'V. (79) 

In the case of invariants up to the fourth degree, M < 4, we see that the closure equations (|66l) and (fHT]) are 
themselves expressed only in terms of the Kahler potential. In these cases, to get a complete solution of the 
problem, it remains only to solve an integrability condition for K.{w,'w). 

In the cases of linear and quadratic invariants, the integrability condition is linear and we have the general 
solution depending respectively on one and two arbitrary real functions (see [7]). In the higher degree cases, 
the integrability conditions are nonlinear and a general solution is lacking: moreover, in general, they give 
rise to an overdetermined system for the unknowns A and VE". This is consistent with the fact that only 
isolated cases of integrable Hamiltonian systems with invariant of degree higher than two are known. Failing 
to work in full generality, a useful approach is to make an ansatz for one of the unknown (e.g. A) and solve 
for 4". In this way several new integrable and superintegrable systems have been found with a cubic ([8]) 
and a quartic ([9]) second invariant. 



4.5 The leading order term of strong invariants 

Looking for a strong invariant, so that the integrability condition must be solved for arbitrary values of E, 
it turns out that the function Gm can no longer be arbitrary. In our previous works we have shown its 
polynomial structure with degree equal to M for M < A and guessed that this is the generic behavior for 
any value of M. In this section we prove that this conjecture is true by explicitly computing the function 
with the aid of the already known solution in the Cartesian case. 

In [7] for the linear and quadratic cases we have found respectively the conditions 

^{S[{z)} = 0, (80) 

3{5^'(^)} = 0. (81) 



In [8], for the cubic case, the condition 



^^^Ss}=0 (82) 



was found, whereas in for the quartic case, the condition 



r4 



S^^54^=0 (83) 



was found. Therefore we can guess the general conditions 



M 



d 

3? <! -^Sm }=Q, M odd, (84) 



M 



d 

'^{^Sm\ ^0, M even. (85) 
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Conditions from ([50]) to (|83p have been obtained working on the integrabihty condition for the Kahler 
potential K. Since it is not possible to write an explicit expression for the integrabihty condition at arbitrary 
M, we prove ([H^ and ([55)) by explicitly calculating them. The easiest way to proceed is to exploit the 
solution in Cartesian coordinates and after that performing the transformation to complex coordinates: the 
reason for not working directly in the complex frame is that, in (j55p and (j56p . part of information on the 
structure of Sm is hidden in the lower order coefficients through the energy constraint: this information could 
be recovered only having the solution of the complete system of equations for Ck up to the integrabihty for 
Co. 

Comparing ([55)1 and ([55]) with (|55|) we get the relation 

M 

SM{x + iy) = Cm=Y. i^^"'"^(m,Af)(a;,y). (86) 

As it is natural, the leading order complex coefficient depends only on leading order Cartesian coefficients. 
Considering l|54p with A: = M + 1, we get the equation for the leading order Cartesian coefficients, namely 

dxA^,-i,M)+dyA^3M)^^, j = 0,l,...,M+l, (87) 
whose solution can be immediately found (see Hietarinta [5], sect. 3.0): 



A(ra,M) = Ty^a(j+k,r,%-kM)^ y 

k=Q j=n 



where a(^j^k,m-k,M) are real integration constants. 
From the relation 

d.^^{d,-idy) (89) 
we have the expression for the r — th derivative 

dl^^Y , A -ifdl-'dl (90) 

In order to compute the operator ([90]) we need the following intermediate result 

ar^a^^"/ = ^^J^^ ,),(fc!',), ^'''^^'^^''^ for A; > s and J > r - s (91) 



whereas the result is zero if fc < s or j < r — s. Using it, the action of the operator ([90| on the monomial 
x^y^ is 

-jr|:Mr (,.:l),,,(,_,/!^;(,_.), ^^-"v--. (92) 

Actually, we are interested only in the in the action of d^^: in this case, since each derivative corresponds to 
decreasing the degree of the monomial by one, only the highest degree terms, with j -\- k = M , survive the 
action of . Applying ([92]) we get: 



QM^M-u^k l_ V( M\ k\{M-k)\ 

y ~ 2^^ {M-s)\s\{s-k)\{k-s)\'' y ■ ^^^^ 



s=0 



Analogously to those in (PT|) , we have that both conditions s — fc > and k~ s>Q must be satisfied. They 
implies just s = A; so that (j93|) turns out to be simply 

dTx^'-'y'^ = ^i~^fjM^^km-k)\ ^ {^^f^M\. (94) 
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The same result holds exchanging k with M — k: 

d^'xV'-' = (95) 
Let us now denote with ^v/) the highest degree part of A(^„^^]\^y. 

A / I 

Aim,M) = (-1) (M-m)!m! ^ ' ^ ^ 

Using (|94p we then have 



M! 



2 



Remembering (j86p we finally get 



M 



dfSM = ^ **^-'"af A(„,M) (98) 

m=0 

JH, ,-ni /\/r|2 

M 

= ^M!a(M,o,M) E (M-m)!m! (^^^^ 

m=0 ^ ' 

which is just what we wanted to prove since it is equivalent to (j84p and (|85p . 

5 Examples 

We illustrate the applications of the general approach described so far with a selection of results, many of 
which are new. In particular, we provide some weakly integrable systems with linear and quadratic invariants 
and a recipe for higher order examples. In order to compactify formulas, from hereinafter with the subscript 
we denote the partial derivative with respect to the corresponding variable (except when used to denote a 
component of momentum). 

5.1 Weakly integrable systems with linear invariants 

The ansatz is 

Ii=Sp + Sp, (102) 

where for simplicity we have suppressed the subscript in the S function. The system of equations ensuing 
from the conservation condition is the following: 

= 0, (103) 

5R{(G5)J = 0. (104) 

Equation p03p agrees with confirming that S can be an arbitrary analytic function, 

S = S{z). (105) 

According to (j69p . the confornial transformation is given by 

dz 

— ^F'{w) = S{z{w)) (106) 
aw 
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or equivalently 

. = x + .y^/^. (107) 

Introducing the "conformal" potential 

G = \F'\'^G = \S\^G = SSG, (108) 

ea. (|104p reduces to 

^{G^} = 0, (109) 

which is readily solved in 

G = g{Y), (110) 

where g is an arbitrary real function and, according to the definition of the coordinate transformation, Y is 
the imaginary part of w. In this coordinates, the invariant assumes the 'normal form' 

h=P + P (111) 

or simply 

h=Px. (112) 
Equation (|104p . in view of (|105p and recalling the definition of G, can be rewritten as 

^{S'{E-V)- SV^} = Q. (113) 

We recognize the structure of eq. (f28|l . where now 

h = 3?{5'(z)} = (114) 

and 

/o = mSV),} - (115) 
and we see that (|114p coincides with ((50)) and solution (jllOp now applies to V, so that we have 

g{Y{x,y)) 

with g arbitrary. The general form of the second invariant in the original real coordinates is 

h=^{S}p,+^{S}Py. (117) 

We may provide two interesting classes of weakly integrable systems admitting linear invariants. The first 
is obtained by the simple observation that, if we choose the level surface E' = 0, it is no longer necessary that 
condition (|114p be satisfied. Any analytic function S = S{z) provides a solution through the corresponding 
conformal transformation. If y, as above, denotes the new coordinate 

then the solution is given by the pair (|116H117"|l . provided we consider only motions at the energy level E = 0. 

The second class of weakly integrable systems is obtained with the following trick. Let us consider the 
analytic function f{z) and consider then the conformal transformation (|106p with S{z) given by 

S= y (119) 

a + f'[z) 
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with a constant. Recalling definition (|108p . let us consider the "flat" conformal potential G = \. In this 
case, relation (|108|) . using (|119p gives 

G = E-V=^ = a' + a{f + /') + |/f . (120) 

We can therefore interpret as the fixed value of the energy constant and get as a consequence the family 
of potentials 

V{z,z;E) = -VEir + f)-\ff. (121) 

To complete the solution, we have to write explicitly the coordinate transformation generated by pi9p . that 
is 

"^"^ =az + f{z). (122) 



S{z) 



The invariant is still of the form (|117p . 



5.2 Weakly integrable systems with quadratic invariants 

The ansatz is 

l2 = Sp"^ + Sp'^ + K. (123) 
The system of equations ensuing from the conservation condition is the following: 

S, = 0, (124) 
Kg + SG, + ^S'G = 0. (125) 

Equation (|124p is already familiar. Since K is real, equation (|125p has the following integrability condition 

^{S"G + 3S'G, + 2SG,,} = 0. (126) 

However, as above, we can directly simplify eq. ()125p . According to ([69|) . the conformal transformation is 
now given by 

— = F'{w) ^ v/^R^ (127) 



dz 



(128) 



G^\F'\'^G = \S\G = ^/s^G, (129) 



or equivalently 

With the conformal potential 
eq. (|125p reduces to 

Kto + G^^ 0, (130) 
that is the first example of the set ([75)1 . Its integrability condition is 

3{G^«,} = (131) 

which is readily solved by 

G = A{w + w) + B{w-w), (132) 
where A and B are arbitrary real functions. Putting this solution in (|125p gives 

K = B -A. (133) 
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Comparing with (|129p and using real "separating" coordinates, 

X^^{w}, r = (134) 

the solution for the original function G is then 

A{X;E) + B{Y;E) 
G= ^ . (135) 

The invariant takes the 'normal form' 

I2 = P^ + P^ + B-A=^{P^-P^) + B-A (136) 
or, in the original real coordinates, 

h = ^Re{S)ipJ ~ py^) + lm{S)p,py + B{Y{x,y)) - A{X{x,y)). (137) 
Equation ()131|) can be rewritten as 

^{S"{E - F) - 3S"K - = 0- (138) 

We again recognize the structure of ea. ipS)) . where now 

A = ^{S"{z)} - (139) 

and 

/o = Q{S"V + 3S"14 + 2SV;:^} = (140) 

and we see that ()139p coincides with (|81|). It can be proven (see '28', section 2) that this condition, in 
addition to warrant strong integrability, is also necessary and sufficient in order that the conformal factor 
15*1 can be written as a sum of two functions of X and Y, say 

\S\ = Sx{X) + Sy{Y). (141) 

Since from (I140p we have that solution (I135|) now applies to V, we get 

with suitable A and B. In this case, X and Y are properly referred to as separable coordinates. From (|135p 
and (|14ip . we deduce the relations 

A{X;E)^ESx-A{X), (143) 
B(Y;E) ^ ESy - B(Y), (144) 



so that, eliminating the energy parameter through the Hamiltonian 

Sx + Sy 



E^H^^^^l±^n±A±£, (145) 



the general form of the strong quadratic invariant is 



J _ Sy{PI+2A)-Sx{P^+2B) 
Sx + Sy 



(146) 
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We mention here some interesting examples of weak integrability with quadratic invariants. Consider 
first the polynomial function S{z) — iz^ . This is the simplest polynomial which gives a potential which is 
not automatically integrable at arbitrary energy. The corresponding conformal transformation is given by 

w = 2-1/2(1 -i)lnz (147) 

or in terms of the real variables using polar coordinates 

X = i=(0 + lnr) (148) 

y = ^(e'-lnr). (149) 
v2 

From the relation (|135p it then follows that the potential given by 

y = r-2[A(re^) + B(re-'')], (150) 
is integrable at zero energy for arbitrary functions A and B. In [28j, the class of systems with 

A(X) = i(C-sin\/2X), (151) 

= ^(C- sin\/2r), (152) 

where C is a real constant, has been investigated. Using polar coordinates, the explicit form of the potential 
is 

Vir.e)^'^^:^^^^^ (153) 

and that of the second invariant is 

l2iPr,Pe,r,9) ^ rprPe - cos9 sin{\iir) (154) 

where 

rpr = xp^ + ypy (155) 
Pe = xpy-ypx- (156) 

A natural question one can ask is if the system is still integrable at arbitrary energy values. As it is well 
known, it is possible to answer this question by performing careful numerical investigations, but, on purely 
analytical grounds, the problem is in general quite difficult. In [28], this potential has been proven to be 
non-integrable at values of the energy — ei < E < €2, with small enough ei,2, for every value of C in the 
interval (0, 1]. The proof has been obtained by applying the Poincare non-existence theorem of additional 
invariants (see, e.g. [311 [32]), working in the transformed coordinates X,Y, taking as unperturbed part of 
the Hamiltonian the integrable system at i? = and as a small perturbation the term associated to the 
conformal factor. 

A second example is that given by a generating function of the form 

S{z) = z-^-. (157) 

The conformal transformation can then be written in terms of the real variables as 

X ^ ^r^cos2e ^ ^{x^ -y^), (158) 

y = ir2sin26' = xy. (159) 

The corresponding potential is 

V{X,Y) ^ r^[A{X) + BiY)]. (160) 
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since the conformal factor is 



Choosing for example 



1 



\SiX,Y)\^^^j=== = -. (161) 



A{X) = (162) 
B{Y) = {2 + a)Y^ + b, (163) 

where a is a constant such that < a < 2 and b > 0, we get the family of potentials 

V{x,y) = r^{x^ + y* + ax^y^ + b). (164) 

This potential has a relative maximum at the origin where V{0, 0) = 0, absolute minima placed symmetrically 
in the four quadrants around the origin and grows as for r —^ oo. It allows bound motion for all the 
admissible values of the energy. For a = 2 it is rotationally symmetric and therefore it is "superintegrable" 
at zero energy. The second invariant is 



h{Px,Py,x,y) = ^ {pI~pI) " ^PxPy + x'^ + y'^ " {4 + a)x^y^ . (165) 



What is remarkable in this case is that already at energies slightly below or above the integrable level E = 0, 
the dynamics is strongly chaotic. This shows that identifying weakly integrable systems does not necessarily 
provide a nearly integrable system. 

5.3 Weakly integrable systems with higher-order invariants 

Exploiting the results obtained in ^ i9j , several families of weakly integrable systems admitting cubic and 
quartic invariants can be obtained. Both cases can be described with a similar approach. The invariants 
have the form 

h = 2Re{S3P^ + R3p}, (166) 
h = 2Re{ S'4P'* + i?4P^ } + Co . (167) 

The conformal transformations (|69p respectively give the normal forms 

l3 = 2Re{P^ + R3P}, (168) 
h^2Re{P^ + RiP^} + Co, (169) 

with 

^ 2-Af 

i?Af = V Rm, M = 3,4, (170) 
in agreement with definition ([72|l . Introducing the Kahlcr potential, eq. ([77|) provides the solutions for Ri\i: 



Rm = -^IC^w- (171) 
In the cubic case, the closure equation (j66p expressed in terms of the Kahler potential gives 

?Stmwn,)C^M = 0. (172) 



In the quartic case, the closure equation is that determined by eq. (|75p that, with j = 2 and the notations 
adopted here, gives 

{Coh = -^{R4&h. (173) 
2G 
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The reality of Co requires the integrabiHty condition 



Equations p72ll74p iUustrate the important difference between the previous hnear and quadratic cases 
{M — 1,2) and the higher order {M > 2) cases. In the linear (see ea. (|109l) ') and the quadratic case (see 
ea. (|13ip l. the integrabihty conditions are linear differential equations: their solution is in terms of one or 
two arbitrary functions respectively. Now, equations (|172p and (|174p are both nonlinear partial differential 
equations: they appear very complicated and in practice impossible to solve in full generality. Only isolated 
systems can be identified. Using the general expression ([75)l . we get in both M = 3 and AI — 4 cases a 
quadratic equation of the form 

hE^ + hE + /o = 0, (175) 

where the fk, k — 0,1,2 depend on S, A, 5" and their derivatives. Looking for strong invariants, S is 
determined as above (cfr. eg. ([82)) and ea. l[83|) ) and several solutions can be obtained by suitable assumptions 
on A and ^' ^8, ^ . 

Looking for weak invariants, a large class of solutions can be obtained in the simplest eventuality E = 0. 
In fact, in this case the only condition we have to satisfy is /o = and it can be proven that this equation 
is a nonlinear PDE in the variable 4* only. In fact, we get respectively 

5R{(*«,^*«,^)^} = 0, M = 3 (176) 

and 

5{(*u;«,«;«'«,tD + 2*.„^.ffi^'^„),„} = 0, Af = 4. (177) 

Now, any solution of these equations in whatever coordinate system can be used to generate a new weakly 
integrable system (at E = 0) using an arbitrary conformal transformation not included in the families 
and IMl). 



6 Conclusions 

Much of the structure and ideas presented in this work has come from the geometric formulation using the 
Jacobi metric as described in [7] and further developed in [8] and [9] . One of the important results of those 
earlier works was the unification of the notions of weak and strong integrability in a common framework. 
Another was the identification of the crucial role of the conformal transformations for analyzing integrability 
in two dimensions. However, while the geometric picture has been very useful in these and other respects, 
it is clear from the present work that an approach which is more closely related to the more traditional 
Hamiltonian picture can also be very effective. 

The conformal transformations have provided the clue for finding explicit conditions for the leading 
order term of strong invariants of arbitrary degree. The conditions, as given in equations ([84]) and ([85]) . are 
natural generalizations of the conditions for low order invariants. The examples given in section 5 provide 
an illustration of the apparent ease by which it is possible to construct interesting weak invariants. We hope 
that this work can serve as a starting point for improving the understanding of higher order invariants of 
both the strong and the weak type. 
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